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We demonstrate dynamic stabilization of an unstable strongly interacting quantum many-body 
system by periodic manipulation of the phase of the collective states. The experiment employs a 
spin-1 atomic Bose condensate initialized to an unstable (hyperbolic) fixed point of the spin-nematic 
phase space, where subsequent free evolution gives rise to squeezing and quantum spin mixing. To 
stabilize the system, periodic microwave pulses are applied that manipulate the spin-nematic many- 
body fluctuations and limit their growth. The range of pulse periods and phase shifts for which the 
condensate can be stabilized is measured and the resulting stability diagram compares well with a 
linear stability analysis of the problem. 



It is well-known that unstable equilibria of physical 
systems can be dynamically stabilized by external peri- 
odic forcing [1]. The inverted pendulum stabilized by 
vertically vibrating the pivot point ('Kapitza's pendu- 
lum') provides a classic example of this non- intuitive phe- 
nomenon and was first demonstrated over 100 years ago 
[2, 3]. Dynamical stabilization has a broad range of ap- 
plications including rf Paul traps and mass spectrometers 
[4], particle synchrotrons [5], and optical resonators [6]. 
Standard analysis of the pendulum problem [7] begins 
with the linearized equation of motion y = (a + /3(t))y, 
where y is the angular deviation from the vertical orien- 
tation, a > 0, and /3(t) is the periodic forcing term. For 
= 0, the linearized equation corresponds to an inverted 
harmonic oscillator, and the perfectly inverted position 
y = y = is an unstable (hyperbolic) fixed point charac- 
terized by exponentially diverging and converging orbits 
in phase space. For /3 ^ 0, determining the stability 
condition is a straight-forward problem of linear stability 
analysis. The case of harmonic forcing, f5(t) = (3cos(uit) 
at frequency u> yields the well-known Mathieu equation 
[8], which among other things describes the motion of 
charged particles in Paul traps. 

Recent advances in ultracold atomic physics provide 
opportunities to investigate unstable equilibrium phe- 
nomena of interacting quantum many-body systems fea- 
turing well-characterized and controllable Hamiltonians 
[9]. By changing the dimensionality of the system, tun- 
ing the interaction strength [10], or magnetically quench- 
ing a spin system [11], it is possible to study excita- 
tions across a quantum phase transition described by 
generalizations of the Kibble-Zurek mechanism and ap- 
parent relaxation to non-thermal steady states in a well- 
isolated quantum system [12]. Beyond these fundamental 
issues, non-equilibrium dynamics can generate squeezed 
states and non-Gaussian states that are potential re- 
sources for quantum enhanced measurements [13] and 
quantum information processing [14]. Dynamic stabiliza- 
tion of non-equilibrium many-body Bose-Einstein con- 
densates (BEC) has been suggested by tuning the sign 
of the interaction strength [15-17] and by time-varying 
the trapping potential in a double well BEC [18, 19]. 
Here, we demonstrate dynamic stabilization of the non- 
equilibrium dynamics of a multi-component spinor con- 



densate. The spin dynamics following a magnetic quench 
are well-described by an unstable inverted harmonic os- 
cillator in the spin-nematic degrees of freedom [20] . The 
non-equilibrium dynamics are stabilized by periodic ap- 
plication of phase shifts to the collective states of the 
system. Viewed broadly, the experiment demonstrates 
open-loop discrete control of a many-body quantum sys- 
tem, and is conceptually related to spin decoupling or re- 
focusing techniques used in nuclear magnetic resonance 
(NMR) [21], bang-bang control of non-interacting two- 
level quantum systems (qubits) in quantum information 
processing [22] and the quantum Zcno effect [23-29] . This 
work paves the way for control of non-equilibrium many- 
body systems via direct manipulations of the quantum 
fluctuations. 

Experimental concept 

The experiment uses small spin-1 rubidium-87 conden- 
sates containing just a single domain, such that the dy- 
namic evolution occurs only in the internal spin degrees 
of freedom. Within this limit, the spinor order parameter 
of the condensate can be written as ip = (C+i, Co, C-i) T 
where Q = ^fple l9i , and p. L = \Q\ 2 are the fractional pop- 
ulations of the three Zeeman states of the / = 1 87 Rb 
ground state hyperfine manifold. The symmetry of the 
spin and magnetic field interactions conserves the con- 
densate magnetization m = {p\ — p-i) and constrains 
the populations p±i = (1 — po ± m)/2. For the m = 
case that is relevant for these experiments, the spin dy- 
namics arc determined by the Hamiltonian [30]: 

7i= ^[.T 2 -(1-.T 2 )C0S^] + |(1-X) (1) 

with x = 2p a — 1 and 6 S = + i + #_i — 29 a . The strength 
of the spin interaction and quadratic Zeeman energy are 
characterized by c and q, respectively, and c < for fer- 
romagnetic 87 Rb. This Hamiltonian has the form of a 
classical non-rigid pendulum and is similar to that for 
the double- well Bose-Hubbard model [30, 31]. The con- 
densate has an unstable equilibrium point at x = 1 for 
small quadratic Zeeman energies and undergoes a quan- 
tum phase transition at q = 2\c\. The nature of the 
instability is more clearly revealed in the spin-nematic 
phase space of the Cartesian components of the spin 
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FIG. 1: Illustration of the dynamic stabilization method, (a) The condensate is initialized at the pole of the spin- 
nematic Bloch sphere, x — 1,S± = Q± = 0. The condensate has Heisenberg-limited uncertainties in S± and Q±. (b) Initial 
free evolution of the condensate produces spin-nematic squeezing along the diverging manifold of the separatrix. (c) The 
quantum state of the condensate is quickly rotated to the converging manifold of the separatrix using a microwave pulse, (d) 
Subsequent free evolution unsqueezes the condensate, returning it close to the original state, (e) Continued free evolution again 
generates spin-nematic squeezing, (f) Long term stabilization is achieved by repeating the (c,d,e) sequence (blue line) whereas 
the unstabilized condensate rapidly evolves away (red line). 



vector Si and nematic (quadrupole) tensor Qij [32]. In 
this space, the spin dynamics are represented on a Bloch 
sphere with axes {S±,Q±,x} where Sj_ = S 2 + S 2 . and 
Q\ = Q%z + Qyz- This sphere is shown in Fig. 1(a) 
together with the dynamical orbits of the system for 
q < 2|c|. The unstable fixed point is located at the pole 
of the sphere at the intersection of the two manifolds of 
the separatrix that divide the space into phase-winding 
and oscillatory phase orbits. In the neighborhood of the 
pole (x ~ 1), the Hamiltonian can be written: 



2 (g-2|c|) 2 

!±- H A -L' 



(2) 



with corresponding canonical equations of motion: Sj_ = 
qQ± and Qj_ = (2|c| — q) Sx- These can be combined 
into Q± — (2\c\q — q 2 ) Q±, which, for q < 2|c|, describes 
an inverted harmonic oscillator that is inherently unsta- 
ble. 

The experiment uses a condensate initialized with the 
atoms in the / = l,m/ = state that corresponds to a 
coherent state centered at the pole of the Bloch sphere 
(x = 1,S± = Q± = 0) as shown in Fig. 1 (a). Evolu- 
tion from this state generates Gaussian squeezed states 
at early times (localized near the pole) [32] and a rich 
variety of non-Gaussian states at later times as the sys- 
tem evolves on the full Bloch sphere [20]. The evolution 



from the initial state is driven by quantum fluctuations, 
ASj_ = AQj^ = TV -1 / 2 , where N is the number of atoms 
in the condensate. These fluctuations are depicted in the 
figure as a Husimi-Q distribution. Initial evolution in the 
neighborhood of the hyperbolic fixed point reduces the 
fluctuations along the converging manifold of the separa- 
trix and grows the fluctuations along the diverging man- 
ifold of the separatrix; together, these create squeezing 
in the spin-nematic phase space. 

The continued growth of the quantum correlations 
eventually destabilize the system and lead to evolution 
away from the fixed point [20]. Dynamic stabilization 
can be achieved by preventing the build-up of these cor- 
relations. In our experiment, this is accomplished using 
periodic phase shifts of the spinor order parameter that 
manifest as a rotation in the S± , Q ± plane about the ori- 
gin. The concept is illustrated in Fig. 1 (a-d). During 
the initial evolution to a squeezed state, the squeezed 
uncertainty ellipse aligns along the diverging manifold of 
the separatrix. In order to prevent further growth of the 
correlations, the state is quickly rotated to the converg- 
ing separatrix using a microwave pulse. Further evolution 
'unsqueezes' the condensate and returns it approximately 
to the initial state, thereby stabilizing the system. Long- 
term stabilization is realized by periodic repetition of this 
cycle. Although we discuss the technique using a rotation 
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FIG. 2: Experimental demonstration of dynamic sta- 
bilization of spin dynamics, (a) Free (unstabilized) evolu- 
tion of the spin population, po, due to quantum spin mixing 
(po = (x + l)/2). (b)-(e) Stabilization of the spin dynamics 
to the po = 1 unstable equilibrium using periodic rotations 
of the spin-nematic quadrature phase with the periods indi- 
cated. In each case, microwave pulses are used to rotate the 
quadrature phase through an amount A 6* = — jir rad. The 
open circles correspond to the mean measured value and error 
bars indicate the measured standard deviation. The red line 
and blue shaded region are smoothed interpolations to guide 
the eye. Note the reduced vertical scale in (b) and (c). 



angle corresponding to the angles between the manifolds 
of the separatrix, the condensate can be stabilized for a 



range of angles as will be shown. 

Observation of dynamic stabilization 

The experiment begins with a condensate containing 
N = 32,000 atoms initialized in the \f = 1,771/ = 0) 
hyperfine state held in a high magnetic field (2 G). To ini- 
tiate spin dynamics, the condensate is rapidly quenched 
below the quantum critical point by lowering the mag- 
netic field to 210 mG. In order to stabilize the dynamics, 
the {S±,Q±} quadrature phase is periodically rotated 
by an angle A9 with a period r. The rotation is imple- 
mented using 2ir pulses on the |/ = l,m/ = 0) <-> |/ = 
2, to/ = 0) clock transition that effectively shift the phase 
of the |/ = 1, to 



/ = 0) spinor component by an amount 
A8 = 7T (1 + A/Vl + A 2 ) where A = S/fl is the detun- 
ing normalized to the on-resonance Rabi rate [32]. This 
pulse effects a quadrature rotation A8 — —A8 . Finally, 
the spin populations of the condensate are measured by 
imaging the condensate following expansion in a Stern- 
Gerlach magnetic field gradient to separate the to/ spin 
components. The atoms are probed for 100 /is by three 
pairs of orthogonal laser beams, and the fluorescence sig- 
nal is collected by a CCD camera. 

The experimental results demonstrating dynamic sta- 
bilization of the condensate are shown in Fig. 2 for sev- 
eral different pulse periods and are compared with unsta- 
bilized free evolution of the spin dynamics. The evolu- 
tion of the unstabilized condensate exhibits a long initial 
pause (~100 ms) in the evolution of (po) (during which 
spin-nematic squeezing is generated) followed by large 
amplitude oscillations that are characteristic of a pendu- 
lum initialized close to the unstable point. The results 
for the stabilization are shown for pulse periods ranging 
from 30-60 ms; in each case, the microwave pulse in- 
duces a quadrature rotation of — |tt rad in the S±,Q± 
plane that is close to the angle between manifolds of the 
separatrix. The data show that the stabilization is close 
to perfect for pulse periods <40 ms for evolution times 
out to 0.5 s. For longer pulse periods, the stabilization is 
still effective for some time, but then eventually degrades. 

In Fig. 3, the measurement of the evolution of the 
transverse spin fluctuations, AS± is shown for different 
pulse periods. To measure S±, a resonant RF pulse is ap- 
plied to the / = 1 manifold to rotate S± into S z before 
the trap is turned off. With no stabilization, the fluctua- 
tions grow exponentially, because Sj_ is not aligned to the 
minimum quadrature squeezing angle, which is approxi- 
mately the converging manifold of the separatrix. When 
the condensate is stabilized, the calculations show that 
fluctuations of S± first increase then decrease during each 
pulse cycle, which reflects the squeezing and unsqueezing 
of the condensate. The data show the expected periodic 
evolution of the fluctuations and also show a significant 
reduction of the fluctuations compared with the unstabi- 
lized condensate. However, the overall level of the mea- 
sured fluctuations are higher than predicted, likely due 
to atom loss induced by the microwave pulses. 

In Fig. 4, we show that the condensate spin dynam- 
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FIG. 3: Transverse magnetization. The fluctuations of 
the transverse magnetization, Sj_, are measured (solid lines) 
for both free evolution of the condensate and for the dynami- 
cally stabilized dynamics for two different pulse periods. The 
results are compared to a theoretical calculations (dashed 
lines). The dB line corresponds to the N^ 1 ^ 2 standard 
quantum limit. 

ics can be stabilized for time scales comparable and even 
longer than the 1/e lifetime of the condensate. The life- 
time of the stabilized condensate is ~500 ms, which is 
reduced from the 1.2 s lifetime measured without the 
microwave pulses. The shorter lifetime is possibly due 
to imperfections in the microwave Rabi pulses used to 
phase shift the condensate that leave some population in 
the / = 2 state, which is short-lived in the trap. The 
stabilization works quite robustly, despite this additional 
loss. We have performed two additional checks to verify 
that the loss does not significantly affect the dynamics of 
the system. In the first, we have studied the evolution of 
the condensate under periodic pulses with A8 = n (the 
periodicity of the phase space) and verified that the con- 
densate undergoes normal spin mixing. In the second, 
we have turned the stabilization pulses off after 500 ms 
and verified that the system again undergoes normal spin 
mixing. 

Determining the stability region 

We now turn to an investigation of the stability diagram 
for the condensate. Although it is conceptually simplest 
to understand the stabilization in terms of periodic evo- 
lutions along manifolds of the separatrix of the phase 
space, the condensate can be stabilized with a range of 
phase shifts and periods (A9,t). The range of stability 
is measured for quadrature phase rotations AO € [— it, 0] 
and microwave pulse period r from 5 ms to 100 ms. 
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FIG. 4: Long time stabilization dynamics of the BEC 
condensate. The fractional population, po, (red, left) is 
compared to the total population of the condensate (blue, 
right). The evolution of the condensate is stabilized for peri- 
ods much longer than the condensate lifetime of ~ 500 ms. 



For each combination of (A#,t), the spin population is 
measured after 195 ms of evolution, where the unstabi- 
lized condensate shows the maximum spin-mixing (see 
Fig. 2). The results of the measurement are shown in 
Fig. 5, where each measurement point is the average of 
three experimental runs. For shorter period pulses, the 
condensate is stabilized with a wide range of quadrature 
phase shifts. For long period microwave pulses, the range 
of quadrature phase shift capable of stabilizing the dy- 
namics shrinks and reaches an asymptotic value close to 
the angle between the separatrices, AO — cos _1 (— 1 — -). 

The results are compared with a theoretical stability 
analysis shown in solid lines. Defining the state vec- 
tor of the system y = {Sj_,Q±) T , the evolution of a 
full stabilization cycle can be written in matrix form as 
M = R[A6>] • exp[mr]. R is the matrix for a rotation 
about the origin in the S±, Q± plane, and A9 = A0 s /2 
is the amount of the quadrature phase shift induced by 
the microwave pulse, m is the matrix form of the equa- 
tions of motion (y = my), and r is the period between 
pulses. The stability condition such that (S±,Q±) stays 
bounded to the origin is |Tr[M]| < 2, which is shown as 
solid lines in the figure. The green lines show the sta- 
bility envelope using the measured values for the spinor 
dynamical rate c = —6.5(5) x 2tt Hz and the magnetic 
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FIG. 5: Stability diagram. To map the stability region, 
the mean po is measured at the time of the maximum spin 
mixing for the unstabilized condensate. The solid curves are 
the envelope of the calculated stability region using linear 
stability analysis, while the dashed curves use a time-averaged 
Hamiltonian approach. The green curves use the measured 
values B — 210 mG, c = —6.5 x 2tt Hz, and the red curves 
use -B e ff = 150 mG, c = —6.5 x 2n Hz. 

field B — 210(10) mG that determines the quadratic Zee- 
man effect q — 71.6 x B 2 Hz/G 2 . Although the measured 
stability region is in good qualitative agreement with the 
theoretical envelope, there is an offset in the asymptotic 
phase for long periods between experimental and theory 
which is about 0.3 rad. This discrepancy is larger than 
can be accounted for in the uncertainties of the deter- 
mination of c and q, the former determined by direct 
measurement of coherent spinor oscillations of the con- 
densate, and the latter determined by RF and microwave 



spectroscopy on magnetic field sensitive hyperfine transi- 
tions. Using a magnetic field value of B c g — 150 mG for 
the theoretical calculation provides a much better agree- 
ment with the data, as is indicated by the red solid lines. 
We suspect that small differential AC stark shifts of the 
Zeeman states induced by the 850 nm optical trap beam 
(at the level of a few Hz) are responsible for this discrep- 
ancy. The measurements also seem to indicate that the 
system exhibits stability outside the analytical envelope 
for small phase shifts, A9 > — 7r/4. This is an artifact of 
the fact that the measurement is determined at a rela- 
tively early time (195 ms), where, for small applied phase 
shifts, the evolution of the condensate is similar to the 
freely evolving case with smaller amplitude spin-mixing 
oscillations and large non- Gaussian fluctuations. 

It is also possible to derive a stability diagram using 
a time-averaged effective Hamiltonian method. For c < 
0, the mf = state is stable if the quadratic Zeeman 
energy is greater than 2|c| or less than zero. For the 
applied phase shifts using the periodic microwave pulses 
we can calculate an effective quadratic Zeeman energy, 
q c s — q + hA9/r. Due to the cyclic nature of the phase 
space, the quadrature shift pulse can be interpreted as a 
positive or negative pulse. Because of this ambiguity of 
the phase shift direction, q c g is double valued everywhere. 
For both q e g > 2\c\ and q c g < there is no longer a 
hyperbolic fixed point centered on the pure mf = state, 
and hence wherever these conditions are met for both of 
the q c g values, the system should be stable. The stability 
boundaries determined using this approach are shown by 
the dashed lines in Fig. 5. They are more restrictive 
than the linear stability analysis, although they agree 
asymptotically for shorter pulse periods. 

Conclusion and outlook 

These experiments demonstrate new methods of manip- 
ulating out-of-equilibrium quantum many-body systems, 
drawing together ideas from classical Hamiltonian dy- 
namics and quantum squeezing of collective states. In 
future work, it should be possible to engineer highly 
entangled states of the system (including Schrodinger 
cat-like states) using extensions of these ideas [33] and 
to demonstrate other fundamental phenomena such as 
Berry phases [34]. 
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In this Supplementary Information, we derive the sta- 
bility regions shown in the main paper for pulsed quadra- 
ture phase shifts applied periodically during the evolu- 
tion. The derivation follows standard methods of sta- 
bility analysis for step-wise evolution used for example 
in optical resonator stability analysis. We compare this 
method with an analysis based on the time-average of 
the Hamiltonian that provides more restrictive condition 
in disagreement with the experiment. 



I. SEMI-CLASSICAL EQUATIONS 

Within the single mode approximation (in which the 
spin states share a common spatial wavefunction) , the 
mean-field order parameter is represented by a complex 
vector (Ci,CojC-i) T where £i represents the amplitude 
and phase of the classical field for the mode associated 
with the 7Ji f state given in the index. A mean field analy- 
sis of the spin Hamiltonian [1-3] produces the dynamical 
equations [4]: 



reduce to [5]: 
2c 
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po and 9 S are canonically conjugate variables and the 
dynamical phase space is defined by the corresponding 
spin energy functional of the condensate: 



£ = + c Po 



(1 - Po) + Vi 1 - Po) 2 -m 2 cosO, 
+pm + q(l - po). 



(2) 



m and 9 m are also canonically conjugate variables. How- 
ever, 9 m is cyclic because it does not appear in the energy 
functional, and thus m is conserved. 



II. DYNAMICAL ANALYSIS 



iKi = EiCi + c[(pi +po- p-i)Ci + CoC-i] 

ih(o = E ( + c[(pi + p_i)Co + 2CiC-iCo] 

= S-iC-i+c[(p-i+Po-Pi)C-i + CoCi] 

where c is the spin interaction energy, pi = Ni/N = \Ci\ 2 
is the fractional population of the rrif — i component, 
and Ei is the magnetic field energy for each mode. 

A convenient parameterization of the order parameter 
is given by 
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where 6 m = 0\ — 9—± is the magnetization or Larmor 
precession phase, 9 S = 9\ +9-\ — 29o is the spinor phase, 
and m = (Ni — N_i)/N = p 1 — p_ 1 is the fractional 
magnetization, which is a constant of the motion. 

Using this parameterization, the dynamical equations 



The first step is to recast the energy functional Eqn. 
2 in terms of the spin and nemantic operators. We will 
restrict ourselves to the m — case which is relevant for 
the experiments. In this case, Eqn. 2 becomes: 

£ = cp (1 - po) (1 + cos0 s ) + q(l - po). 
Substituting x = 2p — 1 gives 

£ = C -(l-x 2 )(l + co S 9 s ) + ^(l-x) 

which is the same as Eqn. 1 of the main paper up to 
constants. 

Next, we evaluate the expectation values of the several 
of the mean field spin and quadrupole operators (pro- 
vided in Appendix) using the parameterization for the 
order parameter in Eqn. 1 with m = 0, 



(S x ) = 2y / p (l - po) cos y cos -y 

(Qyz) = 2Vpo(1 - po)shiy cos-^ 

(S y ) = -2v/p (l -p )cosy sin^y 

(Qxz) = -2i/po(i -p )siny siny 
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Note that for the mean field all of these quadrature val- 
ues range from —1 to 1. From these we calculate the 
relationships 



Si 



= S 2 x + S 2 y =Ap Q (l-p )cos 2 ^ 
= 2p o (l-po)(l+cos0 s ) 



(i 



COS I 



.) 



A 



= Qlz + Q 2 yz = 4po (1 - Po) sin 2 
= 2p (1 - Po) (1 - cos0 s ) 



1 



1 - x 2 ) (1 - cos0 s ) 



Si + Ql + x 2 = 1. 

Substituting these relationships back into the energy 
functional gives 

For low depeletion x — > 1 corresponding to the top of the 
Bloch sphere, we can write the following approximation: 



£ = 



'I, 



si. 



(3) 



For the / = 1 87 Rb condensates, c < 0, and so for 
large positive values of q, the system is isomorphic to 
a stable harmonic oscillator, while for small q, the sys- 
tem is isomorphic to an unstable ('inverted') harmonic 
oscillator. The cross-over point at q = 2\c\ is a quantum 
phase transition. 



III. STABILITY CONDITION 

From Eqn. 3, the canonical equations of motion are 
derived as follows 



S± 



2 d£ 
frdQ± z 
2 dS 



hdS ± 



qQ±z 

---(2c + q)S ± , 



where c = c/Ti and q — q/% are angular frequencies. This 
can be written in a matrix form 



S± 



which is more amenable to the stability analysis to fol- 
low. Defining this matrix as m, the full dynamics from 
one pulse to another including the quadrature phase shift 
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FIG. 1: Stability mapping. Simulation values of po at the 
time of the maximum spin mixing are used to map the sta- 
bility region and analytical fit 



(which is just a plane rotation in S± and Q± z ) is given 
by 



M = R[A0] • exp[rm] 



where r is the period between pulses, R is a rotation ma- 
trix, and A9 = A9 s /2 is the amount of the quadrature 
phase shift. From this point the stability condition is ex- 
actly the same as an optical resonator using ray tracing 
techniques. The determinant of the one pass evolution 
matrix must meet the condition |Tr[M]| < 2. This in- 
equality is used to mark the boundaries of the analytic 
stability region, which compares favorably to simulations. 



IV. EFFECTIVE QUADRATIC ZEEMAN 
ANALYSIS 



Another way to analyze the stability is through the 
use of a time-averaged, effective quadratic Zeeman. For 
c < 0, it is known that the system does not evolve from 
the mf = state if the quadratic Zeeman energy is 
greater than 2|c| or less than zero. We note that the 
quadratic Zeeman portion of m is just a two dimensional 
rotation matrix rotating the quadratures by an angle of 
qr/h. The microwave pulse generates an instantaneous 
quadrature rotation of A9. Therefore over one period of 
the periodic microwave pulse sequence we can calculate 
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an effective quadratic Zeeman, q e n = q + HA9/t. The 
phase space is cyclic with a periodicity of 7r in quadra- 
ture phase. So an instantaneous phase shift of A9 is 
equivalent to a phase shift of AO — ir, which is a negative 
contribution to q c g. Because of this ambiguity of the 
phase shift direction, q e g is double valued everywhere. 
Both values are needed for the time-averaged analysis of 



the stability. For both q c g > 2\c\ and q e g < there is no 
longer a hyperbolic fixed point centered on the m,f = 
state. Wherever these conditions are met for both of 
the q e g values, then the system is trivially stable. This 
region is bounded by the dashed lines in the stability dia- 
gram. The dynamical solutions asymptote to these lines 
for short time between pulses. 



Appendix: Spin and Quadrupole Operators 



TABLE I: Spin-1 dipole operators. Expectation values of these operators are components of the angular momentum vector. 
Matrices in spherical polar basis |/, m/). 



Sy 



s z = 




S x = (a\a + aja-i + oj&i + a^aoj 
S y = a\a Q — ojo-i + ajai + a^ao) 

S z — (&\ai — a^d-ij 



TABLE II: The spin-1 quadrapole operators. Expectation values of these operators are moments of the symmetric traceless 
quadrapole tensor. Matrices in spherical polar basis \f,mf). 



Qyz 




Qxy 
Qxx 



Qyz = ^75 (-a\a + aja_i + ajai - a^ao) 
Qxz = ^75 (a\ao — aja_i + al&i — al^o) 
(-afa_i + a^a^ 



Qyy = — jja^a+i + fajao — loljO-i — a^O-i — a^xa+i 
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